A072  792 


QUEEN  MARY  COLL  LONDON  (ENGLAND)  DEPT  OF  MATERIALS  F/G  11/9 

FRACTURE  MECHANICS  OF  INELASTIC  POLYMERS. (U) 

APR  79  EH  ANDREWS  DA-ERO-75-G-OA9 


1 


UNCLASSIFIED 


FACULTY  OF  ENGINEERING 
QUEEN  MARY  COLLEGE 


«. -I..-—  --  • ..  • ■ .r-i.. 


FRACTURE  MECHANICS  OF  INELASTIC  POLYMERS 


Final  Technical  Report 
by 

E.H.  ANDREWS 
April  1979 


EUROPEAN  RESEARCH  OFFICE 

United  States  Army 
London  England 

GRANT  NUMBER  DA-ERO  75-G-049 

Department  of  Materials 
Queen  Mary  College,  London 


Approved  for  Public  Release;  distribution  unlimited 


79 


FRACTURE  MECHANICS  OF  INELASTIC  POLYMERS 


| E.H.  Andrews 

I 

Department  of  Materials,  Queen  Mary  College,  London  El  4NS , UK. 


: 


i 


Abstract  - The  author's  Generalized  Theory  of  Fracture  Mechanics 
has  been  applied  to  several  quite  different  polymeric  materials 
which  have  in  common  a non-linear  and  inelastic  deformation 
behaviour.  Thus  the  normal  Linear  Elastic  Fracture  Mechanics 
cannot  be  applied.  When  applied  to  experimental  data  the 
theory  provides  a unique  characterisation  of  fracture  resistance 
and  also  sheds  light  on  the  molecular  or  microscopic  fracture 
mechanisms  of  the  material  concerned . 


INTRODUCTION 


Fracture  mechanics  (FM)  is  concerned  with  the  propagation  of  existing  cracks 
and  sets  out  to  define  the  conditions  for  such  propagation.  The  first  FM 
theory  was  that  of  Griffth  [ 1]  who  proposed  a criterion  for  crack  growth 
(and  thus,  by  implication,  fracture  of  the  body)  based  upon  energy 
conservation.  He  said  that  a pre-existing  crack  of  length  c in  a stressed 
body  would  propagate  if 
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Here  the  LHS  is  the  'energy  release  rate',  being  the  decrease  in  the 
elastic  stored  energy  (£)  of  the  system,  with  unit  area  of  growth  in  the 
crack  surface  area  (A) , and  y is  the  surface  free  energy  of  the  solid. 
Using  linear  elasticity  theory  he  was  able  to  evalulate  the  energy  release 
rate  for  a plate  containing  a short  crack  under  a uniform  tensile  stress 


a . Thus,  for  plane  stress, 
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where  E is  the  Young's  modulus.  Equating  the  RHS  of  equation  (2)  to  y, 
he  obtained  the  critical,  or  fracture, stress  as, 
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Following  Griffith,  other  workers  recognised  that  materials  were  not 
generally  elastic  in  the  high  stress  regime  around  the  crack  tip,  and 
suggested  that  y be  replaced  by  y , a term  including  enerqy  losses 
dissipated  by  inelastic  deformations  in  the  bulk  of  the  specimen  but 
notionally  assigned  to  unit  area  of  crack  surface.  This  same  approach  also 
characterises  the  work  of  Rivlin  and  Thomas  [2]  who  used  T the  'tearing 
energy’  instead  of  2y  . 


An  alternative  fracture  mechanics  was  developed  by  Irwin  [ 3]  who 
characterised  the  intensity  of  the  stress  field  around  a crack  by 
stress-intensity  factor  K 
is  given  by 


such  that  the  stress  0_  at  any  point  P 
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where  f . . are  known  functions. 
ID 


(4) 


Irwin's  approach  is  based  on  linear  elasticity  and  is  commonly  referred  to 
as  "linear  elastic  fracture  mechanics"  (LEFM) . We  shall  use  the  initials 
LEFM  to  embrace  any  FM  theory,  including  Griffith's,  which  relies  on 
linear  elastic  behaviour  in  the  solid. 


In  LEFM,  a small  degree  of  inelastic  behaviour  can  be  accommodated  by 
e.g.  treating  a small  plastic  zone  as  an  extension  of  the  crack.  Large 
inelastic  deformations  cannot  be  treated  by  LEFM,  however. 

Since  most  structural  materials  are  selected  for  their  toughness,  they 
are  often  neither  linear  nor  elastic.  Furthermore  they  violate  the 
linear  elasticity  requirement  that  strains  be  infinitesimal.  For 
rubberlike  solids,  Rivlin  and  Thomas  [2]  derived  equations  for  the  energy 
release  rate  or  tearing  energy  T for  a variety  of  specimen  geometries 
without  appeal  to  linear  elasticity.  Their 's  is  essentially  a non-linear 
elastic  analysis.  Finally  Rice  [4]  discovered  a path-independent  contour 
integral,  denoted  J,  which  can  be  used  to  characterise  the  energy  release 
rate  in  non-linear  elastic  systems. 

Each  of  these  theories  is  characterized  by  a fracture  parameter  which 
achieves  a critical  value  co-incident  with  the  onset  of  crack  propagation. 
In  Griffith's  theory  for  example  this  parameter  is  the  energy  release 
rate  and  its  critical  value  is  the  solid  surface  energy  y.  Thus,  to 


we  have:- 
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[N.B.  The  factor  h appears  because  A is  defined  as  the  area 
of  crack  surface  which  is  twice  that  of  crack  plane] 


Theories  1-3  are  based  on  linear  elasticity  theory,  whilst  4-5  are 
non-linear  elastic  theories.  None  of  the  theories  consider  inelastic 
materials  since  none  consider  what  happens  on  unloading  the  stress  field 
as  the  crack  propagates.  They  .ire  therefore  all  concerned  with 
initiation  events  and  can  only  be  used  to  describe  propagation  a-  long  as 
the  mat.  rial  is  regarded  as  el.i  t it  . Finally,  only  in  Griffith.'  theory 
is  the  critical  value  parameter  explicitly  related  to  the  physical 
properties  of  the  solid. 


Turning  this  around  we  could  say  that  existing  theories  suffer  from  the 
following  deficiencies. 

i)  Many  assume  linear  elastic  behaviour  not  found  in  tough  solids,  or, 
usually,  in  polymeric  materials, 

ii)  All  assume  elastic  behaviour  (i.e.  no  energy  losses), 

iii)  Consequently  none  address  themselves  to  crack  propagation  which 

involves  unloading  of  some  regions  of  the  stress  field  and  inelastic 
losses . 

iv)  In  all  theories  except  Griffith's,  the  critical  parameter  which 
characterises  fracture  resistance,  is  purely  empirical  and  not 
related  in  any  explicit  way  to  the  physical  or  continuum  properties 
of  the  solid. 

The  author's  "generalized  theory  of  fracture  mechanics"  was  introduced  to 
overcome  these  various  problems.  In  particular  it:- 

i)  Copes  with  non-linear  and  inelastic  materials, 

ii)  Concentrates  upon  the  problems  of  crack  propagation,  which 
necessarily  involve  unloading  and  mechanical  loss . 

iii)  Relates  the  critical  fracture  parameter  explicitly  to  the  continuum 
properties  of  the  solid. 

iv)  Because  of  its  realism  in  these  matters,  provides  equations  which  are 
of  predictive  value.  For  example,  they  predict  the  critical  energy 
release  rate  in  terms  of  continuum  and  physical  properties  and  also 
predict  such  phenomena  as  fatigue  and  brittle-ductile  transitions. 

THE  GENERALIZED  THEORY 

The  theory  has  been  published  elsewhere  [ 5-6] , so  that  we  shall  simply 
summarize  its  results  here.  These  results  reside  basically  in  two 
equations . 

Andrews ' first  equation 

The  first  equation  actually  represents  a family  of  equations  giving  the 
apparent  energy  release  rate  in  terms  of  the  applied  constraints  and 
dimensions  of  the  specimen.  These  equations  take  the  general  form 

- d£/dA  = k,  (c  ,[  b /£] ) £ w (5) 

ion  o 

lere  the  LHS  is  the  apparent  energy  release  rate,  i.e.  the  energy  release 
rate  that  would  apply  if  the  material  were  perfectly  elastic.  It  is  thus 
identifiable  with  *5  Q for  linear  solids  or  for  non-linear  elastic 
materials.  On  the  RHS,  1 is  a linear  dimension  and  £ , W are 
respectively  the  strain  tensor  and  the  input  energy  density  in  some 
identifiable  region  of  the  specimen.  Usually  this  region  is  chosen  as 
one  where  £ , W are  uniform  e.g.  regions  remote  from  the  crack  in  Simply 
stressed  specimens.  The  function  k depends  on  c but  also  upon  all 
non-infinite  dimensions  b of  the  specimen  other  ?han  _£ , normalised  by  the 
linear  dimension  £_. 

Typically,  £ is  the  crack  length  if  this  is  short  compared  with  the  major 

dimensions  of  the  specimen.  For  an  edge  crack  of  length  c_,  for  example, 

and  b » c 
n — 

- df/dA  = k.  (C  ) c W (6) 

1 o o 

However,  for  specimens  with  crack  lengths  c >>  b where  b,  is  one  of  the 
major  dimensions  of  the  specimen,  £ would  normally  be  chosen  as  b . 
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Andrews 1 second  equation 

The  second  equation  of  the  generalized  theory  is  an  equation  for  the 
critical  apparent  energy  release  rate, 

-(df/dA)  E 7 = 7 $(c,  T,  E ) (7) 

c o o 

Here  the  actual  value  of  -df/dA  at  any  phase  of  crack  propagation  is 
denoted  by  the  symbol  y , and  is  equated  to  the  'surface  energy'  y 
multiplied  by  a 'loss  function'  $. 

The  'surface  energy'  y is  more  strictly  defined  as  the  minimum  energy 
required  to  break  unit  area  of  interatomic  bonds  across  the  fracture 
plane.  Its  precise  physical  significance  is  investigated  further  in 
what  follows . 

The  loss  function  $ is  dependent  on  the  overall  state  of  strain  of  the 
body  (characterized  by  e ) and  also  upon  any  variable  of  the  system  which 
affects  the  mechanical  losses  of  the  material.  The  most  obvious  of 
these  variables  are  the  temperature  and  strain-rate.  The  latter  being 
controlled  largely  by  the  crack  velocity  c. 

The  explicit  form  of  <*>  is 

$ = k.  (e  )/  {k,  (e  ) - £ 3g  6v}  (8) 

1 ° 1 ° £r 

where 

*1  = ('  1^  ) fg6V 

p is  a fraction  dependent  on  the  curvature  of  the  stress-strain 
curve  of  the  material.  For  linear  solids  p = 0.5. 

g is  a distribution  function  of  the  energy  density  W throughout 
the  specimen . 

S V is  the  element  of  reduced  volume,  namely  the  real  volume 
element  divided  by  H3 

3 is  the  hysteresis  ratio  (i.e.  the  fractional  energy  loss  in  a 
stress  cycle  at  a point  in  the  stress  field) . For  elastic 
solids  0=0  and  for  perfectly  plastic  solids  3=1. 

P denotes  summation  over  the  entire  stress  field. 

PU  denotes  summation  only  over  points  which  unload  as  the  crack 
propagates  i.e.  where  dW/dc  is  negative. 

The  purpose  of  this  report  is  to  detail  the  progress  that  has  been  made 
in  our  understanding  of  the  fracture  of  polymeric  solids  by  application 
of  the  generalized  theory.  In  particular  we  wish  to  emphasize  its  value 
in  characterising  highly  non-linear  and  inelastic  materials  and,  secondly, 
in  relating  their  fracture  properties  to  their  molecular  constitution . 

CHARACTERISATION  OF  NON-LINEAR  AND  INELASTIC  POLYMERS 

The  value  of  LEFM  for  brittle  solids  lies  almost  exclusively  in  its 
ability  to  specify  a single-valued  'fracture  toughness'  for  a material 
independent  of  geometry  and  loading  conditions.  For  inelastic  solids, 
characterisation  has  been  attempted  in  terms  of  J and  , but  with  only 
partial  success.  The  problem  lies  in  the  fact  that  neither  nor  6 


have  explicit  physical  significance.  When,  therefore,  variations  are 
observed  in  these  quantities  it  is  impossible  to  say  if  this  is  due  to  a 
failure  of  the  criterion  or  to  some  genuine  change  in  the  mechanism 
of  crack  growth. 

Generalized  FM  corrects  this  situation  by  defining  a 'fracture  energy'  or 
critical  apparent  energy  release  rate,  analogous  to  the  £ of  LEFM.  To 
determine  this  quantity,  3 , requires  the  measurement  of  a critical 
loading  condition,  W , a dimension  (usually  the  crack  length)  and 

the  determination  of°tfie  runction  k^(eQ).  The  latter  function  is 
conveniently  obtained  by  methods  previously  described  [6]  involving 
load-deflection  measurements  on  specimens  containing  edge  cracks  of 
differing  length,  c.  Referring  to  Fig.  1, -A£  is  given  by  OAB,  OBC  etc. 
and  the  corresponding  AA  values  by  c^,  C2-C]/  etc.  times  the  sheet 
thickness.  Then,  from, 

-d£'dA  = k,  (C  ) c W (6) 

1 o o 

a plot  of  - — (At/AA)  against  W is  a unique  function  having  the  slope 
k.  (£  ) or  k^'ftO  (it  is  more  convenient  to  use  W^_  as  the  variable  since 
ex,  R are  in  any  case  uniquely  related) . The  value  of  W is  given  by 
area  8ae  t specimen  volume. 

Data  for  some  highly  non-linear  materials  is  shown  in  Figs.  2 and  3.  In 
Fig.  2 are  the  load-deflection  curves  for  a blown  branched  polyethylene 
film  (ICI  Ltd.  XHF  film),  and  in  Fig. 3 the  k^ (W  ) function  for  this  film, 
for  a second  PE  film  (ICI  Ltd.  XDG  film)  and  for  a polyethylene 
terephthalate  film.  The  shape  of  the  k.  (W  ) is  characteristic  for 
yielding  solids,  having  been  observed  also  for  highly  ductile  metals  [6], 
the  peak  in  k,  occurring  at  the  specimen  yield  point. 

When  x is  combined  with  the  observed  critical  values  for  W , and  the 
corresponding  crack  lengths,  a plot  of  k^  against  reciprocal 

crack  length  should  give  a straight  line  through  the  origin  having  a 
constant  slope  J . Such  plots  are  shown  in  Fig. 4 for  the  PE  XHF  and 
the  PET  films.  In  this  diagram  the  reciprocal  crack  length  axis  has 
been  corrected  by  a finite-width  term,  namely  the  reciprocal  of  the 
specimen  wiuth,  b 1 . The  values  obtained  are, 

PE  XiiF  4 .5  kJ/m2 

PET  5.5  kJ/m2 

We  therefore  conclude  that  Generalized  FM  is  capable  of  producing  a 
geometry-independent  fracture  parameter  for  highly  non-lir.ear  and 
inelastic  materials  which  is  completely  analogous  to  the  of  LEFM. 

PREDICTION  OF  FRACTURE  ENERGIES 

It  is  an  obvious  goal  of  FM  to  be  ai  le  to  predict  critical  fracture 
parameters  such  as  J from  a knowledge  of  the  continuum  and  mo.ecular 
properties  of  solids.  It  is  onualiy  important  to  understand  the  physical 
ignificance  of  these  par.imeters  in  terms  of  such  material  properties. 
These  objectives  have  been  brought  very  much  closer  by  the  second 
equation  of  Generalized  FM. 


load/kgf 


W /kg  cm 
o ^ 


/ 


Fig. 4.  Determination  of  3 for 

polyethylene  film  (circles) 
and  polyethylene  terephthalate 
film  (crosses) 
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Referring  back  to  equation  8,  we  have 


k.  (e  ) 

j , y { i — SL_ } 

° k (e  ) - 3g  <5v 

1 O “ 

PU 

Clearly,  if  $ and  £ can  be  evaluated,  we  have  a prediction  of  the 
value  of  jf . A complete  theoretical  solution  is  still  beyond  our  reach, 
unfortunately.  This  is  because  the  energy  density  distribution  function 
g is  not  generally  known  for  non-linear  deformations.  However,  g can  be 
determined  experimentally  from  the  strain  field  around  a crack,  and  3,  the 
hysteresis  ratio  is  also  experimentally  accessible  for  the  general  inelastic 
case.  For  the  special  case  of  an  elastic-plastic  solid  it  is  sufficient 
to  set  3 - 1 in  the  plastic  zone  and  3 = 0 in  the  unyielding  regions  of 
the  specimen. 

The  calculation  of  J depends  upon  the  system  under  study.  For 

elastomeric  materials,  with  a molecular  network  structure,  , can  be 

o 

calculated  by  the  theory  of  Lake  and  Thomas  (see  below) . For  plastically 
deforming  solids  jf  is  expected  to  be  the  energy  to  form  and  fracture 
unit  area  of  craze,  since  any  normal  evaluation  of  21  from  a 
macroscopic  strain  field  would  not  contain  the  microscopical  craze  zone . 
Below  we  give  examples  of  both  elastomeric  solids,  including  soft  polymers 
like  low  density  PE  and  plasticized  PVC,  and  of  a polymeric  glass, 
polycarbonate . 


Elastomers,  polyethylene  and  plasticized  PVC 

Materials  used  were  SBR  and  EPDM  elastomers,  low-density  polyethylene  and 
plasticized  PVC.  Details  of  these  materials  are  given  elsewhere  [ 7]  . 

We  have,  from  equation  8,  with  suitable  abbreviation. 


y/y0  = v(ki  - zj 


VS  =i 


Now  <<  y normally,  so  that  approximately 

InE  = lnk^  - ( 3Q/y  ) 


This  approximation  will  be  very  accurate  until  increases  to  about 

0.3. 

Equation  9 now  allows  us  to  check  the  predictions  against  experimental 
values  of  tf as  will  be  seen  below,  given  a knowledge  of  E,  k^  and  . 

The  parameter  3 was  evaluated  by  the  method  of  Lake  and  Thomas  [ 8]  using 
the  parameters  given  in  Table  1.  Those  authors  give 

If  = (2/37D*5  y £ U Nn /z  (10) 
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TABLE  1 


Parameters 

SBR 

EPDM 

PI 

PVC 

lx  1010(ni) 

5.15 

3.32 

2.50 

2.50 

m X 10“  (N) 

12.6 

6.79 

4.57 

10.2 

7 

1.54 

1.83 

1.83 

2.77 

U x 10*’  (J> 

12.8 

8.78 

6.60 

6.60 

0 (K) 

296 

296 

296 

296 

T) 

1.15 

1.15 

- 

p x 10-“(Nm-3) 

0.901 

0.871 

0.922 

0.987 

MX  10”  (N) 

32.7 

16.3 

- 

7.35 

v X 10'”  (irrJ) 

1.93 

1.93 

C,  x 10(MNm'J! 

1 2.0 

2.0 

1.5  i 0.5 

5.5  t 3.5 

I 

1 


r 


i 


TABLE  XX 

Theoretical  and  Measured  Values  o£  the  Parameter  (Em  ) 


Material 

Theoretical 

| 

Expt (this  work) 

, 

Expt (Lake&Lindley) j 

NR 

28 

- 

40  1 

r ■ 

EBR 

34 

65  + 5 

60 

EPDM 

36 

65  + 5 

p-PVC 

29  + 10 

loo  + 50 

■ 

PE 

63  + 11 

200  + lOO 

_ 

Notes:-  (a)  Experimental  values  always  50%  to  100%  higher 
than  theoretical,  both  for  our  own  and 
independent  data. 

(b)  Large  error  bands  occur  for  non-crosslinked 
materials  PE  and  PVC  but  even  these  data  are 
consistent  with  note  (a) . 
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where  y is  a flexibility  parameter,  £ the  length  of  a monomer  unit,  U 
the  bond  rupture  energy  for  a monomer  unit,  N the  number  of  network  chains 
in  unit  volume  and  n the  number  of  monomer  units  between  cross-links.  Also 

N = 2v  = 2C^/kO  (11) 

where  k is  Boltzmann's  constant,  the  first  i-looney-Rivlin  constant  and  0 
the  absolute  temperature.  The  cross-link  density  V is  the  sum  of  chemical 
links  (v*)  and  physical  entanglements  (v1)  and  its  effective  value  takes 
account  of  chain  ends  thus 

v = (v*  + V1)  (1  - r|p/MV*)  (12) 

where  g is  a constant,  p the  density  and  M the  average  molecular  mass  of 
the  molecules.  In  the  case  of  PE  and  PVC,  V*  = 0 and  all  the  effective 
"cross-links"  are  physical  in  origin. 

The  results  for  CT0  are  shown  in  Table  II  together  with  data  of  Lake  and 
Lindley  [ 9] . These  authors  found  that  their  theoretical  values  were 
consistently  a factor  of  about  two  smaller  than  experimental  values 
obtained  from  -atigue  experiments  and  ascribed  the  difference  to 
deficiencies  in  the  theory. 

The  evaluation  of  £ was  carried  out  as  follows.  For  the  case  of  an  edge 
crack  in  a semi-infinite  sheet,  this  term  takes  the  form 

4 2LT  6 {x  — • + y — } 6x  6y  (13) 

ox  3y 

PU 

where  8(x,y)  is  tne  hysteresis  ratio  (equal  to  the  fractional  energy  loss 
in  a stress-strain  cycle) , x and  y are  the  reduced  variables  X/c  and  Y/c 
respectively  and 

f = W(x,y) /Wq  (14) 

The  summation  is  taken  only  over  points  P which  are  unloading  (PU) 
i.e.  at  which  the  energy  density  W decreases  if  the  crack  propagation. 

in  order  to  evaluate  £ v/e  must  clearly  know  8 and  W as  functions  of  x,y. 
Since  also  a propagating  crack  represents  a dynamic  situation,  and  6,  W 
are  dependent  upon  strain  r _es  as  well  as  strain  levels,  we  need  to  know 
strain  e and  strain  rate  t as  functions  of  x,y. 

To  obtain  these  data  is  laborious,  but  not  intrinsically  difficult.  The 
experiments  carried  out  to  this  end  are  summarized  here. 

a)  An  X,Y  grid  was  printed  on  the  sheet  containing  the  edge  crack  with  a 
mesh  size  0.2mm  (Plate  1). 

b)  The  crack  was  caused  to  propagate  by  deforming  the  sheet  in  an  Instron 
testing  machine,  the  rate  of  crack  growth  being  controlled  by  the 
overall  strain  applied  to  the  specimen. 

c)  The  strain  field  surrounding  the  propagating  crack  was  recorded  by 
photographing  the  grid. 

d)  The  crack  velocity  was  measured  optically. 

e)  Hysteresis  ratio  8 was  obtained  for  the  material  as  a function  of  strain 
level  and  strain  rate  by  independent  simple-extension  strain-cycling 
tests  on  the  Instron. 


Plate  1 Propagating  crack  (depth  1mm) 


in  a sheet  of  EPDM  rubber, 
showing  strain  field 


Typical  hysteresis  loops  are  shown  in  Figure  5.  Ideally,  of  course, 
hysteresis  data  is  needed  in  bi-axial  strain  states,  but  this  could  not  be 
done  simply.  It  is  an  assumption  that  at  a given  temperature  and  strain 
rate  £,  the  value  of  3 is  uniquely  related  to  the  input  energy  density  W 
regardless  of  the  bi-axiality  of  strain. 

f)  The  appropriate  strain  rates  at  a point  P are  calculated  from  the  crack 
velocity  and  the  major  principal  strain  £^  at  P from  the  approximate 
formula 


d£. 


dY 


dY 

dt 


c 


(15) 


where  the  Y axis  is  the  crack  axis  and  c the  crack  velocity. 


At  this  stage  we  have  all  the  data  necessary  for  the  calculation  of  E. 
The  actual  evaluation  is  carried  out  by  the  double  integration  of  the 
argument,  thus 


E 


3 (£,  e)  x dfdY  + 

c 


(e,  i)  YdfdX 


where 


R 

x 


R 

y 


o 


h X dW 

h Y dW 


(16) 


(17) 


To  evaluate  R , plots  are  made  of  energy  density  W as  a function  of  the 
product  hX  for  various  values  of  Y.  W is  derived  from  the  principal 
strains  £ , c7  (obtained  directly  from  the  distorted  grid  taking  into 
account  that  the  grid  lines  are  not  themselves  the  principle  axes  of 
strain),  oy  use  of  the  Moon  ;y-Rivlin  equation. 


In  all  cases  X,  Y are  referred  to  the  undeformed  state. 


R is  then  obtained  by  graphically  integrating  3*  with  respect  to  W up  to 
tiie  peak  in  the  W vs  X characteristic,  since  only  this  region  will  unload 
if  the  crack  propagates.  R is  similarly  evaluated  for  various  values  of 
X and  E obtained  by  graphical  integration  of  the  R vs  Y and  R vs  X 
curves  (see  Figures  6 and  7 which  represent  summation  over  hal^  the 
strain  field  only) . 


It  was 


found  that, 


d X » 


(18) 


the  second  term  being  negative  and  of  the  order  of  5%  of  the  first.  This 
simplified  subsequent  work  which  neglects  the  second  term  and  thus  halves 
the  effort  involved  in  evaluating  E. 


Hysteresis  loops  for  EPDM 
rubber  at  20  C 
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A step  in  the  computation  of 
I (see  text) 
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Fig. 7.  As  Fig. 6 
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Fig.  9 
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Plot  of  Z versus  J 1 for 
polycarbonate.  As  in  Fig. 8., 
Jo  is  the  negative  slope 


Using,  from  equation  8, 


k.  3 

E - *1  - 


we  predict  a linear  relationship  between  1 
slope  of  (-k^  and  an  intercept  k at  3 1 
excellent  agreement  with  this  prediction. 


and  I with  a negative 
•*  0.  Figure  9 shows 


Significance  of  *J_  for  polycarbonate 

The  value  of  derived  from  Fig. 9 is  large,  namely  3.3  kJ/m2 . This 
compares  with  0.2  kJ/m2  for  polyethylene  and  even  less  for  elastomeric 
materials.  It  is  obviously  impossible  to  interpret  'J  for  polycarbonate 
in  terms  of  the  network  theory  used  earlier  to  justify  high  values  in 
elastomeric  and  similar  materials.  Even  if  the  entanglement  network  in 
polycarbonate  is  considered  to  replace  the  vulcanized  network  in 
elastomers,  it  is  difficult  to  see  how  such  a high  figure  could  be  expected. 


However,  when  the  "J  value  found  here  is  compared  to  the  *j  C value 

O r . 

extrapolated  by  Ward  and  co-workers  [ 10J  for  plane  strain  (totally 
brittle)  impact  fracture  in  similar  materials,  a very  close  agreement 
is  observed.  Under  these  conditions,  no  plastic  flow  occurs  but  only 
a microscopical  craze  which  precedes  the  crack.  Since  in  calculating  I 
no  account  was  taken  of  the  very  high  energy  density  gradients  encountered 
in  this  microscopical  craze  zone,  it  follows  that  the  energy  to  form 
the  craze  must  be  included  in  the  value  of  'J  . We  conclude  therefore 
that  the  value  3.3  kJ/m2  obtained  for  ’J  is  ?he  craze  formation  energy. 

SUMMARY 


The  fracture  and  tearing  of  highly  non-linear  and  inelastic  solids, 
ranging  from  elastomers,  through  visco-elastic  film  materials  to 
elastic-plastic  polymer  glasses,  is  well-described  by  the  newly 
introduced  Generalized  Theory  of  Fracture  Mechanics.  Not  only  does  the 
theory  provide  the  ability  to  characterize  fracture  resistance  in 
non-linear,  inelastic  materials,  but  it  goes  further  to  provide  a new 
parameter,  j , which  can  be  related  to  the  mechanism  of  fracture 
i ropagution  it  a molecular  or  microscopical  level.  The  interpretation 
of  0 varies  according  to  the  system  under  investigation.  For  the 
adhesivi  liiu.  e of  a Van  der  Waals  bonded  interface  f ll],  rJ  is  simply 
half  t thi  rinodynumic  work  of  adhesion.  For  cohesive  fracture  of  network 
polymers,  J is  the  energy  needed  per  unit  area  of  surface,  to  break 
single  network  chains.  For  polycarbonate,  'J  is  the  energy  to  form  the 
microscopical  craze.  In  all  cases  the  theory  allows  us  to  relate  the 
macroscopical  fracture  mechanics  to  the  microscopical  physical  processes 
of  energy  loss  (via  8)  and  failure  (via  $ ) respectively. 
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